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ABSTRACT 

00 
O , 

O . One of the most puzzling current experimental physics paradoxes is the arrival on Earth 

^ I of Ultra High Energy Cosmic Rays (UHECRs) with energies above the GZK threshold 

(SxlO^^eV). Photopion production by CMBR photons should reduce the energy of these 
protons below this level. The recent observation of 20TeV photons from Mk 501 (a BL Lac 
object at a distance of ISOMpc) is another somewhat similar paradox. These high energy 
to '• photons should have disappeared due to pair production with IR background photons. A 

common feature of these two paradoxes is that they can both be seen as "threshold anoma- 
lies" : energies corresponding to an expected threshold (pion production or pair creation) are 
^ ■ reached but the threshold is not observed. Several (relatively speculative) models have been 

proposed for the UHECR paradox. No solution has yet been proposed for the TeV-7 para- 
dox. Remarkably, the single drastic assumption of a violation of ordinary Lorentz invariance 
would resolve both paradoxes. We present here a formalism for the systematic description 
of the type of Lorentz-invariance deformation (LID) that could be induced by non-trivial 
short- distance structure of space-time, and we show that this formalism is well suited for 
comparison of experimental data with LID predictions. We use the UHECR and TeV-7 data, 
as well as upper bounds on time-of-fiight differences between photons of different energies, to 
constrain the parameter space of the LID. A model with only two free parameters, an energy 
scale and a dimensionless parameter characterizing the functional dependence on the energy 
scale, is shown to be sufficient to solve both the UHECR and the TeV-7 threshold anomalies 
while satisfying the time-of-fiight bounds. The allowed region of the two-parameter space is 
relatively small, but, remarkably, it fits perfectly the expectations of the quantum-gravity- 
motivated space-time models known to support such deformations of Lorentz invariance: an 
integer value of the dimensionless parameter and a characteristic energy scale constrained 
to a narrow interval in the neighborhood of the Planck scale. 
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1 Introduction 



Significant evidence has accumulated in recent years suggesting that in two different regimes, 
Ultra High Energy Cosmic Rays (UHECRs) and multi-TeV photons, the universe is more 
transparent than what it was expected to be. UHECRs interact with the Cosmic Microwave 
Background Radiation (CMBR) and produce pions. TeV photons interact with the Infra 
Red (IR) photons and produce electron-positron pairs. These interactions should make 
observations of UHECRs with E > SxlO^^eV (the GZK limit) n or of gamma-rays with 
E > 20TeV from distant sources unlikely g, |, |]. Still UHECRs above the GZK limit and 
20TeV photons from Mk 501 are observed. 

Numerous solutions have been proposed for the UHECR paradox (see [^] for a recent 
review). Most of these solutions require new Physics. There are practically no proposals 
concerning the TeV-7 paradox (see however, [B]). It is striking that there are some common 
features in these otherwise apparently unrelated paradoxes. In both cases low energy photons 
interact with high energy particles. The reactions should take place because when Lorentz 
transformed to the CM frame the low energy photon have sufficient energy to overcome 
an intrinsic threshold. In both cases the CM energies are rather modest (~ 100 MeV 
for UHECRs and ~ 1 MeV for the TeV photons) and the physical processes involved are 
extremely well understood and measured in the laboratory. In both cases we observe particles 
above a seemingly robust threshold and the observations can be considered as a "threshold 
anomaly". It is remarkable that in spite of these similarities at present there is only one 
mechanism that could resolve both paradoxes: a mechanism based on the single, however 
drastic, assumption of a violation of ordinary Lorentz invariance. 

The possibility that the cosmic-ray threshold anomaly could be a signal of violation 
of ordinary Lorentz invariance had already been emphasized in Refs. 0, pJP, lp |, Jllll . In 
this work we combine these earlier points with the very recent suggestion ^14 |14|| that 
Lorentz-invariance violation could be the origin of the TeV-7 threshold anomaly, vve analyze 
a general phenomenological framework for the description of the type of Lorentz-invariance 
deformation (LID) that could be induced by non-trivial short-distance structure of space- 
time, and we ask whether there are choices of LID parameters that would simultaneously 
solve the two threshold anomalies while satisfying the constraints imposed by the fact that 
the results of experimental searches [O, |16[ of energy-dependent relative delays between the 
times of arrival of simultaneously emittedphotons are still consistent with ordinary Lorentz 
invariance. We obtain, under these assumptions, strict limits on the possible parameter 
space of LID. The fact that one is at all able to give a quantitative description of both 
threshold anomalies with a simple two-parameter LID model provides encouragement for the 
interpretation of the data as a sign of LID; moreover, it is quite remarkable that the values 
expected from quantum-gravity considerations (most notably the energy scale characterizing 
the deformation being given by the Planck scale) are in agreement with the strict limits we 
derive. 

We review in Sections 2 and 3 the observational background and the theoretical problems 
related to the observations of UHECRs (Section 2) and TeV photons (Section 3). In Section 
4 we describe a special (two parameter) model for LID and we obtain limits on these two 
parameters. In Section 5 we describe a more general five-parameter LID formalism and 
again we constrain the parameter space with the available data. In Section 6 we compare 
our formalism with the Coleman and Glashow ||^ formalism for Planck-scale-independent 
Lorentz invariance violations. We summarize our results in section 7. Appendix A is devoted 
to the K-Minkowski space-time, which is an example of quantum-gravity motivated space- 
time that allows a simple illustration of some of the structures here considered. 



2 UHECRs and the GZK paradox 

The high energy cosmic rays (CR) spectrum depicts a clear break at ~ 5x lO^'^eV. This break 
is accompanied by a transition in the CR composition from nuclei to protons. Above this 
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break the spectrum behaves (with a decreasing statistical certainty due to the small number 
of events) as a single power law N{E) ~ E~^''^ all the way up to 3.2xlO^°eV [0, the highest 
energy CR observed so far. 

A sufficiently energetic CMBR photon, at the tail of the black body thermal distribution, 
is seen in the rest frame of an Ultra High Energy (UHE) proton with E > SxlO^^eV as a 
> 140MeV photon, above the threshold for pion production. UHE protons should loose 
energy due to photopion production and should slow down until their energy is below the 
GZK energy[|. The process stops when CMBR photons energetic enough to produce pions 
are not sufficiently abundant The proton's mean free path in the CMBR decreases 
exponentially with energy (down to a few Mpc) above the GZK limit (~ 5x lO^^eV). Yet more 
than 15 CRs have been observed with nominal energies at or above 10^°±30% eV[|T8l, |19]. 

There are no astrophysical sources capable of accelerating particles to such energies within 
a few tens of Mpc from us (at least not in the direction of the observed UHECRs). Further- 
more if the CRs are produced homogeneously in space and time, we would expect a break 
in the CR spectrum around the GZK threshold: below the threshold we would observe CRs 
from the whole universe; above the threshold we would observe CRs only from the nearest 
few Mpc. The corresponding jump by a factor of ~ 30 — 100 in the extrapolated number 
counts above and below the threshold, is not seen. 

Numerous solutions have been proposed to resolve the GZK paradox (see for a recent 
review). These solutions include, among others, new physics solutions like the decay of 
topological defects, weakly interacting messengers like 5*0 or neutrinos with anomalous cross 
sections at high energies (the 'Z-burst" model). Conventional astrophysics solutions like 
acceleration of UHECRs by GRBs or local AGNs require the ad hoc assumption that Earth 
is located in a not generic place in space-time (we should be nearer than average to a typical 
source by a factor of 5) as well as very strong intergalactic magnetic fields |5D|. Another 



conventional solution, the acceleration of Fe nuclei by magnetars in the galactic neiio, requires 
a new, otherwise unobserved, population of galactic halo objects. Clearly there is no simple 
conservative solution to this puzzle. 

From the point of view of our LID phenomenology it is important to notice that for 
a solution of the GZK paradox it would be necessary (and sufficient) for LID to push the 
threshold energy upwards by a factor of six. In fact, the mean free path of a SxlO^^eV 
proton is almost a Gpc, while the highest observed UHECR energy is 3.2xlO^°eV. 



3 TeV photons from Mk 501 and Mk 421 

HEGRA has detected high-energy photons with a spectrum ranging up to 24 TeV from 
Markarian 501 (Mk 501), a BL Lac object at a redshift of 0.034 (~ 157 Mpc). This obser- 
vation indicates a second paradox of a similar nature. A high energy photon propagating in 
the intergalactic space can interact with an IR background photon and produce an electron- 
positron pair if the CM energy is above 2meC^. The maximal wavelength of an IR photon 
that could create a pair with a 10 TeV photon is 40um. As the cross section for pair cre- 
ation peaks at a center of mass energy of about 3meC , 10 TeV photons are most sensitive to 
30/iM IR photon and the mean free path of these photons depends on the spectrum of the 
IR photons at the ~ 15 — 40;uM range. These wavelengths scale like lOTeV/E for different 
energies. 

^The exact composition of UHECRs is unknown and it is possible that UHECRs are heavy nuclei rather 
than protons. In this case such nuclei would undergo photodisintegration when interacting with CMBR 
photons. The threshold energy for a photodisintegration of a nuclei is several MeV. It just happens to be 
true, purely as a result of a numerical coincidence, that the threshold is reached when the energy of a typical 
nuclei, say Fe, is ^ SxlO^^eV. Thus the GZK paradox is insensitive to the question of what is the exact 
composition of UHECRs. 
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There have been several attempts to model the IR background resulting from different 



cosmological evolutionary models m2, p3l 24, ESl. Recently, new data from DIRBE at 2.2yuM 
Zg, at 60 and 100 /iM [gTf and atl4Dahd 230 /iM and from ISOCOM at 15/iM 



suggest that the IR background is even higher. According to these data the flux of IR photons 
is ~ 2.5xl0~^erg cm~^sec~^sr~^ around 60-120 fiM and falls off by an order of magnitude 
towards 15 fiM. This decrease is important as it would lead to a much shorter mean free 
path for 20TeV photons as compared to the mean free path of lOTeV photons. 

It was originally suggested that the expected break, corresponding to hard-photon dis- 
appearance in the IR background, in the GeV-to-TeV spectrum of AGNs could be used 
to determine the IR background spectrum. This would have been based on searches of a 
distance-dependent break in the spectrum of various AGNs. However, no apparent break is 
seen in the spectrum of MK 501 at ~ 20TeV range, where the optical depth seems to exceed 
unity. Using current IR background estimates Coppi and Aharonian find an optical 
depth of 5 for 20TeV photons from MK 501 (see also [0). This optical depth increases 
rapidly with energy. Thus, photons at these energies are exponentially suppressed, unless 
they somehow evade the pair-production process. 

Unlike the GZK paradox only a few solutions have been proposed for the TeV-7 paradox. 
First, it is possible that there is an upturn in the intrinsic spectrum emitted by Mk 501. 
Such an upturn would compensate for the exponential suppression at this region. Clearly this 
is an extremely fine-tuned solution as the expected energy of this upturn should somehow 
be tuned to the energy at which the optical depth from MK 501 to Earth is unity. This 
energy scale is distance dependent and it puts us in a very special position relative to the 
source. It is of course possible that the IR intensity has been overestimated. A shift in 
the energy estimate of HEGRA would also explain the paradox. Finally Harwit, Protheroe 
and Biermann suggest that multiple TeV photons may be emitted coherently by Mk 501 
and if they arrive at Earth very close in time and space they may be confused with a single 
photon event with higher energy. 

With current data, ~ lOTeV photons from Mk 501 could reach Earth, while ~ 20TeV 
photons are exponentially suppressed. This happens mainly because of the rapid fall off of 
the IR spectrum below 60/im. We conclude that a LID upwards shift of the threshold energy 
by a factor of two would resolve this paradox. 

Having discussed the relevance of Mk 501 for the emergence of the TeV-7 threshold 
anomaly we turn now to TeV photons from Mk 421 (another BL Lac object at a redshift of 
0.031, corresponding to ~ 143Mpc). It is not clear if the spectrum of this source extends 
high enough to pose a paradox comparable to the one indicated by Mk 501. However, we 
note here the simultaneous (within the experimental sensitivity) arrival of ITeV photons and 
2TeV from this source. This was used to limit the time-of-fiight differences between photons 
of different energies to less than 200 seconds. This in turn allowed to establish, through an 
analysis of the type proposed in Ref. [0, an upper limit on Planck-scale-induced LID |1^ 



which will be a key element of our analysis. We call these constraints in the following 
time-of-fiight constraints. 



4 Lorentz-invariance-violating dispersion relation 



We start by considering first, a class of dispersion relations (following 33, 31] for a = 1, 
and ]P3, P3| for a general a) which in the high-energy regime takes thelorm: 
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m, E and p denote the mass, the energy and the (3-component) momentum of the particle, 
Ep is the Planck energy scale {Ep ~ lO^^MeV), while a and t] are free parameters character- 
izing the deviation from ordinary Lorentz invariance (in particular, a specifies how strongly 
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the magnitude of the deformation is suppressed by Ep). Clearly, in (|l|) the "speed-of-light 
constant" c has been set to one. (Note however that in this framework c is to be understood 
as the speed of low-energy massless particles |^.) Also notice that in (|l|) we wrote the 
deformation term in two ways, as a E'^{E/Ep)" correction and as a p'^{E/Ep)°' correction, 
which are equivalent within our present analysis based exclusively on high-energy data, for 
which E p, but would be different when studied with respect to low-energy data. (Of 
course, a given short- distance picture of space-time will have only one dispersion relation; 
for example, in "/t-Minkowski space-time" , the space-time which we describe in Appendix A 
in order to illustrate in an explicit framework some of the structures relevant for our analysis, 
one encounters a deformation of type p'^{E/Ep)°'.) 

In previous works [31-35,9] a slightly different notation had been used to describe this 
same class of deformations, which in particular replaced our t] by two quantities: the scale 
Eqg = |?7|~^/"-Ep and a sign variable ^± = f]/\f]\- The a,T] notation turns out to be more 
suitable for the description of the technical aspects of the analysis discussed here, but it is 
useful to keep in mind that the scale of Lorentz-deformation is obtained as {rjl'^^^Ep. 

As hinted by the presence of the Planck scale, our interest in deformed dispersion relations 
of type (|l]) originates from the fact that such deformations have independently emerged 
in theory work on quantum properties of space-time. We postpone the discussion of this 
motivation to the next Section, where we also clarify which types of generalizations of (|I]) 
could also be motivated by Planck-scale physics. 

While our analysis is motivated by the role that the deformed dispersion relation (|I]) 
might have in quantum gravity, one could of course consider ([^) quite independently of 
quantum gravity^ The quantum-gravity intuition would then be seen as a way to develop 
a theoretical prejudice for plausible values of a and r]. In particular, corrections going 
like (E/Ep)" typically emerge in quantum gravity as leading-order pieces of some more 
complicated analytic structures Ip. This provides, of course, a special motivation 

for the study of the cases a = 1 and a = 2.[f{E/Ep) ^ 1 + ai{E / EpY'^ + ....] Moreover, the 
fact that in quantum gravity the scale Eqg is expected to be somewhere between the GUT 
scale and the Planck scale corresponds to the expectation that rj should not be far from the 
range 1 <ri < 10^°. 



4.1 Deformed thresholds from deformed dispersion relations 

We intend to discuss the implications of Eq. (0) for the evaluation of threshold momenta. 
Before doing that let us briefly summarize the derivation of the equation describing the 
threshold in the ordinary Lorentz- invariant case. Relevant for our phenomenological con- 
siderations is the process in which the head-on collision between a soft photon of energy e 
and momentum q and a high-energy particle of energy Ei and momentum pi leads to the 
production of two particles with energies E2.,E^ and momenta P2-,Pz- At threshold (no energy 
available for transverse momenta), energy conservation and momentum conservation imply 

E, + e = E2 + Es, (2) 

Pi~q=P2+P3; (3) 

moreover, using the ordinary Lorentz-invariant relation between energy and momentum, one 
also has the relations 

2 

q = e, Ei = ^pf + ~ p . + , (4) 
[ 2pi 

■^Having mentioned that of course the deformation (|l|) could be considered independently of its quantum- 
gravity motivation, let us also mention in passing that even outside the quantum-gravity literature there is 
a large amount of work on the theory and phenomenology of violations of Lorentz invariance (see, e.g., the 
37| , which also provide a good starting point for a literature search backward in time). 



recent Refs. R |, 36, 
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where rrii denotes the mass of the particle with momentum pi and the fact that pi (and, as 
a consequence, p2 and ps) is a large momentum has been used to approximate the square 
root. 

The threshold conditions are usually identified by transforming this laboratory-frame 
relations into CM-frame relations and imposing that the CM energy be equal to m2 + m^; 
however, in preparation for the discussion of deformations of Lorentz invariance it is useful 
to work fully in the context of the laboratory frame. There the threshold value pi^th of the 
momentum pi can be identified with the requirement that the solutions for m and p^ as 
a function of pi (with a given value of e) that follow from Eqs. (||), (H) and (P should be 
imaginary for pi < pi^th and should be real for pi > pi^th- This straightforwardly leads to 
the threshold equation 

{m2 + - ml 
Pi,th ^ • (5) 

This standard Lorentz-invariant analysis is modified |^, |ll|, |13], |14|, ^ by the defor- 
mations codified in (|l]). The key point is that Eq. (H) should be replaced by 

l+a ^2 l+a 

e = q + n- , E'^ ~ H + n- — . (6) 

Combining (|^), (0) and (P) one obtains a deformed equation describing the pi-threshold: 

Pu. + nj^^ l^ („, + ,„3),.„ - 1 j • m 

where we have included only the leading corrections (terms suppressed by both the smallness 
of Ep"^ and the smallness of e or m were neglected). 



4.2 Phenomenology 

Early phenomenological interest in the proposal ([Ij) came from studies based on time-of 
flight analyses JO, 
According to (|l[)(ani 
relation v - - ■ - - 



TB| of photons associated with gamma-ray bursts or with Mk 421. 



_ assuming that there is no leading-order deformation of the standard 

dE/dp) one would predict ||3l|, |3^ energy-dependent relative delays between the 
times of arrival of simultaneously emitted massless particles: 



AT 



(a + 1) E" 



p 



where T is the (average) overall time of travel of simultaneously emitted massless particles 
and AT is the relative delay between the times of arrival of two massless particles of energies 
E and E' . The fact that such time delays have not yet been observed allows us to set 
bounds on the a, rj parameter space. In particular, data showing (approximate) simultaneity 
of arrival of TeV photons from Mk 421 were used [|I^ to set the bound |?7| < 3-10^ for a = 1. 



The same data were used in Ref. [0] to set a more general a-dependent bound on rj. 

We combine these existing bounds with the assumption that indeed the UHECR and 
TeV-7 threshold anomalies are due to LID (|lD. The fact that the scale Ep is very high might 
give the erroneous impression that the new term going like pf^h/^p present in Eq. (|^) could 

For given values of 
||35[| (displacement 
old completely 



always be safely neglected, but this is not the case f^. 
a, 1] one finds values of e that are low enough for the ''threshoL 
of the threshold) to be significant 



m |T|, o, 0, |5 



anornaly 

For certain combinations a, rj, e the thres 
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disappears, i.e. Eq. (|^ has no solutions. Assuming (|^) one would predict dramatic depar- 
tures from the ordinary expectations of Lorentz invariance; in particular, if a ~ —7] ~ 1, 
according to one would expect that the Universe be transparent to TeV photons. The 
corresponding result obtainable in the UHECRs context would imply that the GZK cutoff 
could be violated [0 even for much smaller negative values of rj. Positive values of rj would 



shift the thresholds in the opposite direction {e.g. they would imply an even stricter limit 
than the GZK one) and are therefore not consistent with the hypothesis that UHECR and 
TeV-7 threshold anomalies be due to LID (|ID. 

In Figure 1 we provide a quantitative description of the region of the a,-?] parameter 
space which would provide a solution to both the UHECR and TeV-7 threshold anomalies 



while satisfying the time-of-flight constraints |15, M that are still consistent with ordinary 



Lorentz invariance. The curve describing the time-oFnight constraints was obtained using the 
information that there is [0 an upper bound of order 200 seconds to the difference in time of 
arrivals of 2TeV photons and ITeV photons simultaneously emitted by Mk 421 (at redshift of 
0.031). The two threshold-anomaly curves reported in Figure 1 were obtained using Eq. (|^ 
with mi = and m2 = = 5-lO^eV (TeV photons j + 'j e^ + e~ threshold analysis) 
and with mi = m2 = 9.4-10®eV and = 1.4-10'^eV (UHECR p + ^ p + re threshold 
analysis^). In light of the analysis of the experimental situation provided in Sections 2 
and 3, we obtained the UHECRs curve by requiring sufficient LID to explain the factor-6 
threshold shift 5-lO^^eV— 3-10^°eV, while for the TeV photons curve we required a factor-2 
threshold shift lOTeV 20TeV. Even though the shift is more significant in the UHECR 
context, it is the requirement to explain the TeV-7 threshold anomaly that provides a more 
stringent constraint, as one should expect since in our LID, which is motivated by Planck- 
scale physics, the violation of ordinary Lorentz invariance is suppressed by some power of 
the ratio E/Ep. 

Considering the diverse origin and nature of the three relevant classes of experimental 
data that we are considering, the fact that there is a region of the a,-!] parameter space 
consistent with all these constraints is non-trivial, and this in turn provides encouragement 
for the interpretation of the threshold anomalies as manifestations of LID. Moreover, it is 
quite striking that this region of parameter space, in spite of being relatively small, does 
contain one of the two mentioned quantum-gravity-motivated scenarios: a = 1 and 1 < 
7] < 10^. The other quantum-gravity-motivated scenario, the one with a = 2 and 1 < 77 < 
10^, is outside the relevant region of parameter space, being consistent with the absence of 
relative time delays and the UHECR threshold anomaly but being inconsistent with threshold 
anomaly for multi-TeV photons. 

Concerning the consistency of the interpretation of the threshold anomalies as manifesta- 
tions of LID it is also important to observe that the modified dispersion relation (^, in spite 
of affecting so significantly the GZK and TeV-7 thresholds, does not affect significantly the 
processes used for the detection of the relevant high-energy particles. For the significance of 
the threshold modification a key role is played, as evident from equation (|^), by the smallness 
of the energy of the background photons. The effect of (|lj) on atmospheric interactions of 
the relevant high-energy particles is instead suppressed by the fact that in these atmospheric 
interactions the "targets" , nuclei or electrons, have energies much higher than those of the 
background photons. 



■^The dominant contribution to the GZK cutoff actually comes from the A resonance, so one might find 
appropriate to replace the sum of the proton mass and the pion mass with the mass of the A in the UHECR 
threshold formula. However, the difference between niA and nip + to^ would only introduce a relatively 
small correction in our UHECR limit which is not our dominant lower limit (a much stricter limit comes 
from the TeV-7 anomaly). Moreover, once the contribution to GZK from the A is avoided one would still 
have a (weakened) GZK cutoff from non-resonant photopion production and this would anyway lead to the 
limit we use. 
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Figure 1: The region of the a, rj parameter space that provides a solution to both the UHECR 
and TeV-7 threshold anomalies while satisfying the time-of- flight upper bound on LID. Only 
negative values of r] are considered since this is necessary in order to have upward shifts of 
the threshold energies, as required by the present paradoxes. The solid thick line describes 
the timc-of-flight upper bound. The region above this line is excluded. The solid thin line 
and the dotted line describe the lower bound on LID obtained from the present UHECR 
(solid thin hne) and TeV-7 (dotted hne) threshold anomalies. The anomahes disappear in 
the region above the hues. Within the narrow region between the dotted line and the solid 
thick line the time of flight constraint is satisfied and both anomalies are resolved. The two 
vertical segments at a = 1 and at a = 2 (i.e. at l/a = 1/2) correspond to the two favored 
quantum-gravity scenarios. The behaviour of the curves for upper and lower bounds on LID 
with respect to the bottom-left corner of the frame can be understood by noticing that at a 
fixed a ordinary Lorentz invariance can be reached taking the 77 — > limit, while at fixed 77 
this requires taking the a — > 00 (i.e. the l/a — > 0) limit. 
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5 A more general LID formalism 



Having shown that the simple two-parameter family of Lorentz-invariance-violating disper- 
sion relations (|l|) provides a solution of the UHECRs and TeV-7 threshold paradoxes, we 
turn now to a more general five-parameter LID formulation. The motivation for this formu- 
lation comes primarily from theory work on short-distance (so called, "quantum gravity") 
properties of space-time, in which modifications of space-time symmetries are encountered 
quite naturally. In particular, quantum-gravity effects inducing some level of nonlocality 
or noncommutativity would affect even the most basic flat-space continuous symmetries, 
such as Lorentz invariance. This has been recently emphasized in various quantum-gravity 
approaches [31-33,39-50] based on critical or noncritical string theories, noncommutative ge- 
ometry or canonical quantum gravity. While we must be open to the possibility that some 
symmetries are completely lost, it appears plausible that some of them are not really lost but 
rather replaced by a Planck-scale-deformed version. Some mathematical frameworks which 
could consistently describe such deformations have emerged in the mathematical-physics 



literature |B3l p8l p9, 44, 51 . An example of these structures is discussed in Appendix A. 



5.1 The five-parameter formalism 

The fact that a simple two-parameter family of Lorentz-invariance-violating dispersion re- 
lations is consistent with all available data is of course of encouragement for the LID 
hypothesis, but, especially since relevant data are expected to improve rapidly in the coming 
years, it is also important to establish how much room for generalizations of (|1|) is available 
in the general framework of Planck-scale-induced LID. 

One way to generalize ([I|) would involve attributing different independent values of a, rj to 
different particles. We shall not pursue this (however phenomenologically viable) possibility, 
since the focus of the present article is on deformations of Lorentz symmetry which could 
be induced by non-trivial space-time structure, and such deformations would most likely 
treat "democratically" all particles. In any case, it is clear that models attributing different 
independent values of a and 77 to each particle end up having a very large number of free 
parameters and available data will not be very effective in constraining such models. We 
shall come back to this point in Section 6, where we consider the alternative (Planck-scale 
independent) Coleman and Glashow scheme for Lorentz-invariance-violation. In fact, 
that scheme corresponds to the choice a = and an independent value of 77 for each particle. 

Another way to generalize the dispersion relation (Q) is to include other deformation 
terms. In a space-time with some non-trivial structure at distances of order one could 
expect that probes with energy much smaller than Ep should obey a dispersion relation of 
type: 

-p^ -m^ = F{E,p,m;Ep) , (9) 

where F is some general function with units of mass (or energy) squared and such that F ^ 
for Ep —>■ 00. Actually, in studies, such as ours, looking only for the leading correction, one of 
the arguments of F can be suppressed: one makes a subleading error by using E"^ —p^ —m? = 
to express one of the variables in F in terms of the other variables. One could for example 
express F as a function of p and m only: F{p,m; Ep). Moreover, the fact that we are only 
looking for the leading correction in the high-energy regime^ allows us to approximate F 
with its leading (if any) power dependence on Ep and (within a given power dependence on 

''it is perhaps worth emphasizing that the low-energy expansion of F{p, m; Ep) may look quite different 
from its corresponding high-energy expansion. In the high-energy regime {p ^ m) the premium is on the 
leading dependence on p while in the low-energy regime (p <^ m) the leading dependence on m is dominant. 
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Ep) leading dependence on p: F{p,m] Ep) ~ rjp^^'^ ^m^ Ep'^. In the high-energy regime 
there is therefore scope for considering the three-parameter family of dispersion relations 

E^-p'^-m^ ^ r]-p^+'^-"-nf-E-'^ , (10) 

where, of course, it is understood that m°" = 1 whenever cr = 0, even when m = 0. (The 
parameter a has been introduced to characterize the type of dependence of the deformation 
term on the mass m, and therefore in our notation there is the implicit prescription that 
m°" — >■ 1 in the formally ambiguous combined limit a ^ 0, m ^ 0.) 

Besides the structure of the dispersion relation a LID can also affect the law of sum of 
momenta. Since our emphasis is here on the phenomenology of LIDs, rather than on their 
formal/mathematical analysis, we limit our discussion of the motivation for this type of effect 
to the example of non- commutative geometry (the "k- Minkowski space-time" ) considered in 
Appendix A. As that example clarifies, it is natural to consider a two-parameter class of mod- 
ifications of the law of sum of (parallel) momenta Ki + K2 ^ Ki + K2 + 5{KiK2)^^^'^^^'^E~^ . 
For our threshold analyses this corresponds to 

Pi-e^pi-e-d ~pi-e, P2 + Ps ^ P2 + Ps + . (11) 

Ep Ep 

Overall we consider a five-parameter space: a,T],a for the dispersion relation and (3,6 
for the description of possible deformations (0) of the law of addition of momenta. The 
analysis reported in the preceding Section corresponds of course to the a 0, 6 limit 
of this more general five-parameter (a, rj, a, (3, 6) phenomenology. 

As appropriate for the present preliminary status of the experimental situation and the 
fact that the two-parameter phenomenology analyzed in the previous Section turned out to 
give a fully satisfactory description of the data, we shall only provide here a preliminary and 
partial exploration of the enlarged five-parameter space. Our exploration of this parameter 
space will also be more detailed in some directions and less detailed in others. In particular, 
we shall limit our analysis to two classes of scenarios, one with a = and one with cr = 
2,6 = 0. This will be sufficient for a qualitative understanding of how different portions of 
our five-parameter space compare with the present experimental situation. 

Retaining the leading corrections in E~^, the threshold analysis in the general five- 
parameter {a, T], cr, P, 6) LID scenario leads to the threshold equation:0 

^Pi,th ( y/msmg \ / ^ 
2e \m2 + m^) \ Ep ) 

In the following sections we apply this equation to several specific cases. To simplify the 
discussion we provide here explicit expressions for the threshold for photopion production: 

^Note that actually the threshold is not necessarily anomalous; in particular, as we already observed in 
Ref. p5[ |, when a = (3 = l, a = Q and 77 = —5 there is a cancellation and the deformed symmetries lead to 
the same threshold equation obtained with undeformed symmetries. 
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GZK,th 



7x10^9 eV 
e/O.OOleV 



1 + r] 1022-2-10.9a-8.15a [(0.87^+" + 0. 13^+") 1 . 15" - 1 



7xl0i9eV/ 



Ep/10i9GeV 



(13) 



_ r 1Q22.1-8.15/3 EGZK,th \ V ^GZA-,ffe/7xl0^^eV Y 

l7xl0i9eVy 1^ Ep/10i9GeV y 



and for pair creation threshold: 
25 TeV 



E. 



e/O.OleV 



1 + r/ 10 



14.8-7.7o--14.6o 



(^.-a _ .K \ ( E,,tH y-'^ / ig,,,/25 TeV ^ 
V "'°H25TeVy \Ep/10^^GeV ^ 



- ^ ini^-s-i^-^'^ / E^^th y / E^^th/2b TeV V 
V25TeVy 1^ Ep/10i9GeV / 



(14) 



where we found convenient to introduce the "Kronecker delta", here denoted with 5^ to 
differentiate it from our parameter 5, to compactly write this equation consistently with 
our conventions for the mi limit. [In deriving Eq. ([T^ ) from Eq. ([T^) it is necessary 
to take into account that, consistently with the conventions and notations we introduced 
(see, in particular, the comments made immediately after Eq. (0)), in the limit mi 
the term mj must be handled according to m^ ^ if cr 7^ ana according to m^ — 1 if 
(7 = 0. Of course, the reader can verify by direct calculation that this prescription gives the 
correct threshold conditions that follow from Eq. ( p!0D in the two cases (7 = and a 7^ 0, 
and reproduces the threshold condition (^ obtained in the preceding Section (which was 
devoted to the case cr = 0, 5 = 0).] 



5.2 Phenomenology with a = 2, ^ = 

In the case a = 2, 5 = there is no deformed law of addition of momenta and the threshold 
equation takes the form 

p,„. ^ + ( (,„^_^.^^)„_i - rn, j - (15) 

For 0" 7^ the LID term in ([l0| ) vanishes for massless particles. Thus, in general in all 



cr 7^ cases (like the one we discuss in this Subsection) the time of flight constraints JTS], [16 
do not limit the LID parameters. 

The constraints obtainable by interpreting the UHECR and TeV-7 threshold anomalies 
as manifestations of LID suggest that this interpretation is quite unnatural in the case 
cr = 2, (5 = 0. The condition that both threshold be pushed upwards leads to the constraints 
77 > 0, a < 1.195. Moreover, in order to describe the threshold anomaly for multi-TeV 
photons one should also make the awkward requirement 77 > 10^^". Having provided in the 
previous section an elegant solution of the threshold paradoxes using o" = we do not pursue 
further this scenario which appears to require a higher level of fine tuning. 

Scenarios with cr 7^ might regain some interest if there are significant new developments 
in the understanding of the threshold anomalies that will point in this direction. In the 
present experimental and theoretical situation we find appropriate to make in the following 
the assumption that a = 0. 
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5.3 General aspects of the phenomenology with a = 

For a = one is left with a four-parameter space on which significant information can be 
gained by combining data on possible time of flight delays, which will only constrain, through 
the prediction (H), the parameters a, rj, and data on the threshold anomalies, which, through 
(p!2[), are relevant for all four parameters a, t], P, 6. 

It is important to observe that positive (discovery) results on both the thresholds and 
the time delays would allow to determine the values of all four parameters. If eventually 
the mentioned time delays are actually observed, and if they are observed in signals from a 
collection of sources diverse enough to allow the determination of the energy dependence of 
the time delays, we would then be able to use (||) to fix a and t]. Then, knowing a and rj, a 
determination of the thresholds could be used to fix /5 and 6. 

While waiting for these eventual discoveries, one can use the present upper limits on 
LID in relative time delays and (preliminary evidence of) lower limits on LID in threshold 
anomalies to reduce the allowed portion of the four-parameter space. We subdivide the 
discussion of this type of phenomenological analysis in three cases: a < (3, a = P and a > p. 

5.4 Phenomenology with a < (3 (and a = 0) 

The case a < P (and cr = 0) is essentially analogous to the case considered in the preceding 
Section with the two-parameter a, rj phenomenology. In fact, for a < P the threshold 
corrections associated with the deformation of the law of addition of momenta are suppressed 
by factors of order {E/Ep)^~°' with respect to the threshold corrections associated with 
the deformed dispersion relation. The constraints derived for a, 77 in the preceding Section 
would still be valid and, as long as we have only lower or upper limits (rather than definite 
discoveries), no constraint could be put on P,6. 

5.5 Phenomenology with a = (3 (and cr = 0) 

For a = P (and a = 0) the upper limit on time-of-fiight LID still constrains only a, r], but the 
constraints on a, rj obtainable by interpreting the UHECR and TeV-7 threshold anomalies as 
manifestations of LID are weakened by allowing also a deformed law of addition of momenta. 
In practice the parameters a, rj and /?, 6 can in a sense "share the burden" of explaining the 
threshold anomalies. To illustrate this mechanism we show in Figure 2 the constraints on 
r], 5 that are obtained for a = /? = 1 (and cr = 0). 

5.6 Phenomenology with a > (5 (and cr = 0) 

For a > P (and a = 0) the threshold corrections associated with the deformed dispersion 
relation are suppressed by factors of order {E/Ep)°'~^ with respect to the threshold correc- 
tions associated with the deformation of the law of addition of momenta. Therefore the 
interpretation of the UHECR and TeV-7 threshold anomalies as manifestations of LID im- 
poses constraints (lower bounds on LID) on the parameters P, S. As always, the upper limit 
on time-of-fiight LID constrains only a,T]. It is worth noticing that if future data should 
indicate that there is no LID relative time-delay effect but there are LID threshold anomalies 
this scenario with a > P would become favored. 

Figure 3 depicts the limits on P, 6 that follow, when a > P, from interpreting the UHECR 
and TeV-7 threshold anomalies as manifestations of LID. The limits on a, t] due to the upper 
limit on time-of-flight LID are still the same as in Figure 1 (but, as just mentioned, the two 
threshold- anomaly curves in Figure 1 do not apply when a > P). 
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Log(6) 



Figure 2: A two-dimensional slice through the five-dimensional parameter space. Shown 
is the T) < and S > region for a = j3 = 1 and (7 = 0. Analogous considerations (with 
exchange of roles between rj and 5) also apply to the corresponding 77 > and 5 < region. As 
in Figure 1, the thick solid line describes the time-of-flight upper bound, while the tentative 
lower bounds on LID that can be obtained from the present UHECR and TeV-7 threshold 
anomahes are described by the thin sohd fine and the dotted fine respectively. Notice that 
for 5 < 10~^^ the range of allowed 77 values is almost unaffected by S, while values of S such 
that S > 10^ '^ are not consistent with the working hypothesis of the present Article: the 
tentative threshold-anomaly lower bound is higher than the time-of-flight upper bound for 
S > 10^ \ 
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Figure 3: The region of the S parameter space that provides a solution to both the UHECR 
and TcV-7 threshold anomalies for a > f3. Only negative values of 6 are considered since, 
when a > f3, this is necessary in order to have upward shifts of the threshold energies, as 
required by the present paradoxes. As in Figure 1, the tentative lower bounds on LID that 
can be obtained from the present UHECR and TeV-7 threshold anomalies are described by 
the thin solid hne and the dotted line respectively. 
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6 Comparison with the Coleman-Glashow scheme 



Coleman and Glashow [H have recently introduced a different scheme (denoted CG scheme 
hereafter) for violation of Lorentz invariance. Modifying the elementary particles Lagrangian 
they suggest a scheme in which there is a different maximum attainable velocity, c^, for each 
particle. The relevant dispersion relations take the form 

E'-p'cl = m'ct, (16) 

where the index a labels the particle. In the language developed in Sections 4 and 5 these 
dispersion relations (|l^) involve two terms, one with a = and a = and the other with 
a = and a = 2. The particle-dependence of could be described by allowing for a 
different independent value of rj for each fundamental particle. At high energies, in which we 
are interested, the a = 0, a = term dominates and rja = — ^ 2c{c — Ca)- The condition 
a = reflects the fact that the CG scheme is not motivated by Planck-scale physics. The 
possibility for each particle to get its own independent value of rj reflects the fact that this 
scheme is not intended as a description of deformations of Lorentz invariance due to non- 
trivial short-distance space-time structure. (If a deformation of Lorentz symmetry is induced 
by the structure of space-time we expect that it would affect all particles in the same way. 
Such a symmetry deformation might allow for a dependence of the correction terms on the 
mass and the spin of the particle but the parameters of the model should not depend on the 
mass, spin or other quantum numbers of the particles.) 

Using again the language we developed in Sections 4 and 5 one can also give an intuitive 
characterization of the way in which the CG scheme and the scheme considered here are 
alternative to one another as strategies for obtaining threshold anomalies. In fact, in that 
language one could describe undeformed threshold^ as associated with a = 0, cr = 0, 5 = 0, 
independently of the value of rj. This corresponds to the fact that in the CG scheme there are 
of course no threshold anomalies if all c^'s take the same value [ca = c — t]/ (2c)). Threshold 
anomalies are generated in the CG scheme by deforming the threshold conditions in the 
direction that corresponds to keeping a = 0,a = 0,6 = but allowing different independent 
values of Cq for each fundamental particle. On the contrary, in our scheme the threshold 
anomalies are obtained by allowing for deviations from a = 0,a = 0,6 = while keeping a 
single (particle-independent) 77. 

In light of these comments it is not surprising that threshold anomalies within the CG 
scheme take the characteristic "Ca — Cb" dependence. In particular, as already observed in 
Ref . 1^ , the description of the UHECR threshold anomaly requires (together with conditions 
on Ca — Cp) that c-,, — Cp > 10~^^. (c^r and Cp are the Ca's for pions and protons respectively.) 
We observe that a resolution of the TeV-7 threshold anomaly within the CG scheme requires 
the additional condition Ce — c^ > 5-10~^^. This combines with the absence M of vacuum 
Cerenkov radiation by electrons with energies up to SOOGeV in such a way that Cg C'y IS 
bound to 5-10^^® < Ce — Cj < 5-10^^^. There is therefore a relatively narrow range of allowed 
values for Cg — c^ just likeQ we found in Section 4 a relatively narrow allowed region of the 
a,ri parameter space. 

^The fact that there are no UHECR and TeV-7 threshold anomaUes in our scheme for a — 0,a = 0,6 — 
can be easily derived directly from the corresponding dispersion relation. This is also implicit in Figure 1, 
which shows that jryl ^ 00 as a — > 0. Also notice that undeformed thresholds are not only obtained for 
a = 0,(7 = 0,5 ~ 0: the thresholds become undeformed also, for example, in the limit a ^ 00 approached 
keeping a — 0,6 — 0. However the undeformed-threshold point a = 0,a = 0,S = is best suited for a 
comparison between our scheme and the CG scheme. 

''Note however that, while the TeV-7 threshold anomaly is used in both, not all the experimental con- 
straints used by the two phenomenological analysis are the same. In particular, the time-of-flight upper 
bound on LID was not used to establish 5-10^^^ < — < 5-10"^^. 
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One important difference between the two schemes is that in our Planck-scale-motivated 
LID the allowed region of parameter space is found exactly where quantum-gravity intuition 
would have sent us searching for new physics, while in the CG scheme values of Ce — c^, in the 
range 5-10~^^ < Ce — < 5-10~^^ do not have any special significance. Another important 
difference between the two schemes is that while the same a, 77 parameters of our scheme for 
LID are also constrained by UHECR threshold data, in the CG scheme Ce — does not play 
any role in the equation for the UHECR threshold and vice versa. Any future development in 
the UHECR threshold data would leave Cg— unaffected. On the contrary, the plausibility of 
the Planck-scale-motivated LID will be strongly affected by future UHECR threshold data: 
if the lower limit on the threshold continues to be pushed higher the overall consistency 
and appeal of the LID model would increase, while the discovery of the threshold not much 
higher than the present 3-10^°eV lower limit would (unless the TeV-7 threshold anomaly 
is eventually understood as a result of systematic errors) rule out the model considered in 
Section 4. 

While the scheme considered here is more tightly constrained by high-energy data (be- 
cause all high-energy data set constraints to the same few space-time related parameters), 
the CG scheme is constrained more tightly than ours by low-energy data. The parameters 
we considered in the present Article, dealing exclusively with the high-energy regime, are 
practically unconstrained by low-energy data since, as discussed in Section 5, the LID we 
considered here might emerge in quantum gravity as the leading order in the high-energy 
expansion of an analytic function whose low-energy expansion looks quite different. On the 
contrary the CG scheme takes a fixed (energy-independent) value of its parameters Ca and 
therefore high-energy and low-energy data can be combined to obtain stricter limits. 



7 Summary and outlook 

In the present Article we took as working assumption that the UHECR and TeV-7 threshold 
anomalies do not have a simple explanation (whereas, especially for the case of TeV photons, 
it might still be legitimate to explore the possibility that systematic experimental errors 
be responsible for the paradox, and other solutions exist for the GZK paradox) and we 
attempted to test the plausibility of a description of the anomalies in terms of a Planck- 
scale-induced deformation of Lorentz symmetry. The results reported in Section 4 certainly 
indicate that this description is plausible. Had we not been considering such a dramatic 
departure from conventional physics, we would have probably gone as far as stating that the 
LIDs we considered provide a compellingly simple description of the anomalies. We do feel 
that the results of Section 4, also taking into account that there is no other known common 
explanation of the two threshold paradoxes, provide motivation for additional theory work 
on the speculative idea of LID and for additional experimental studies aimed at testing the 
class of Planck-scale-induced LID here considered. 

While presently-available data do not in any way invite one to look beyond the simplest 
two-parameter LID examined in Section 4, in preparation for future studies, especially the 
expected improvement of the experimental input, we have developed in Section 5 a general 
parameterization that may prove useful for future attempts to constrain (even rule out) 
Planck-scale-induced LID. We have emphasized the fact that Planck-scale-induced LID, since 
it should reflect the structure of space-time, can be characterized by a small number of 
parameters. In the high-energy regime we found that a very general description of (the 
leading effects of) Planck-scale-induced LID only requires five parameters and we described 
how the determination of a few thresholds together with measurements of the speed of very- 
high-energy particles could fix all five parameters. While we have considered a very general 
class of LID, it should be stressed that we postponed to future studies the analysis of an 
important class of further generalizations of Planck-scale-induced LID deformation 
terms involving a dependence on polarization/spin of the particles. 
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In closing we would also like to emphasize the fact that the experimental data here 
considered represent an important sign of maturity for the general programme of "Planck- 
length phenomenology" p4| , |55| . Whether or not Planck-scale-induced LID turns out to 
successfully describe future experimental data, the fact that at present we are confronted with 
experimental paradoxes whose solution could plausibly involve the Planck length, and that 
certainly the relevant class of observations will eventually be able to rule out various pictures 
of the short-distance (possibly quantum) structure of space-time, shows that, contrary to 
popular folklore, some experimental guidance can be obtained for the search of theories 
capable of unifying gravitation and quantum mechanics. This confirms the expectations. 



that emerged from the general quantum-gravity studies reported in Refs. ||34 



which were based on analyses of the sensitivity of various classes of experiments [52 



51J,|53 



53, 54, 55 



and from analogous studies, primarily focusing on the hypothesis that the unification of 



gravitation and quantum mechanics should involve non-critical strings, reported in Ref. |56 



We thank Daniele Fargion and Glennys Farrar for many informative discussions on UHECRs. 
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Appendix A: /t- Minkowski space-time 

In order to illustrate in an explicit framework some of the structures relevant for our analysis 
of LID, in this Appendix we give a brief descrition of the "/t-Minkowski" non-commutative 
space-time, which was developed in Refs. [Q, |39|, The simplicity of ^-Minkowski, 

which is basically an ordinary Minkowski space-time on wnich however one postulates that 
the time coordinate does not commute with the space coordinates ([xi,t] ~ Xi/Ep), renders 
it very useful for the purpose of illustrating the new conceptual elements required by space- 
times with a nontrivial short- distance structure. 

A first point that deserves being emphasized is the connection between flat nontrivial 
space-times and quantum gravity. In quantum gravity one has the general intuition [Rl, 
that ordinary classical commutative space-times should emerge from some more fundamental 
underlying picture. To very compact (Planckian-energy) probes space-time should look 
completely different from an ordinary classical space-times. On the contrary probes of very 
low energy should not be affected in any noticeable way by the nontrivial short-distance 
structure of space-time. In the intermediate regime (mid-energy probes []34| , ^ ) one would 



expect to be able to use roughly the same language of ordinary classical space-times, but with 
the necessity to introduce some new concepts (such as the little element of noncommutativity 
of K-Minkowski) reflecting the leading-order effects of quantum-gravity at low energies. This 
hierarchy of regimes is to be expected not only in high-curvature space-times (where classical- 
gravity effects are stronger), but also in space-times that appear to be trivial and flat to very- 
low-energies probes. ^-Minkowski is a model (toy model?) of how a probe of relatively 
high energy could perceive a space-time that instead appears to be trivially Minkowski to 
probes of very low energy. 

In K-Minkowski a relation of type (I]) can be obtained as a direct consequence of the 



K-Poincare invariance 33 



|3^, pUf of this space-time. ^-Minkowski therefore provides 

an example of the mention"ed scenario in which an ordinary symmetry is violated but there 
is no "net loss of symmetries" (the 10-generator Poincare symmetry is replaced by the 10- 
generator k- Poincare symmetry). It is in order to capture the essence of these situations 
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that one introduces the terminology "symmetry deformation" (in alternative to "symmetry 
violation" which could be reserved for cases with a net overall loss of symmetries). In 
Sections 4 and 5 we denominate our scheme as a LID just to emphasize that the equations 
we use do not necessarily reflect a loss of symmetry (whether or not they do imply a net loss 
of symmetry depends on the underlying algebraic structures that lead to those equations in 
a given space-time picture). 

Importantly, consistency with the non-commutative nature of K-Minkowski space-time 
also requires |3^, 

This modification emerges at the level of the K-Poincare (Hopf) algebra, and of course 
requires physical interpretation (particle momenta in a noncommutative space-time are a new 
concept). A prescription suitable for handling the ambiguities due to the non-commutative 
nature of ^-Minkowski space-time was given recently in Ref. [p^], and in the cases here of 
interest, which always involve the sum of parallel momenta oitwo particles (at threshold 
particles are produced at rest in the CM frame), it reduces (in leading order in E~^) to the 
prescription that the sum of momenta Ki and K2 can be handled with ordinary algebraic 
methods upon the replacement Ki + K2 Ki + K2 + 5K1K2/ Ep, where 5 is a parameter 
analogous to 77. In the analyses reported in Sections 4 and 5 this would imply — e — > 
Pi — e — 6pie/Ep and P2 + Ps P2 + Ps + Sp2P3/Ep, and actually, since of course we have 
been here only interested in the leading E~^ effect and e -C p2 ~ Pa ~ Pi, one can neglect 
the term of order pie/Ep while retaining the term of order p2Pz/ Ep. 



44, M that the law of addition of momenta be accordingly modified. 
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